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Abstract 

The Hurst exponent H of long range correlated series can be estimated by means of the De- 
trending Moving Average (DMA) method. A computational tool defined within the algorithm is 
the generalized variance cr 2 DMA = l/(N - n)^Jy(i) - y n (i)] 2 , with y n (i) = l/nYj k y{i - k) the 
moving average, n the moving average window and the dimension of the stochastic series y{i). 
This ability relies on the property of a 2 DMA to scale as n 2H . Here, we analytically show that o 2 DMA 
is equivalent to Cnn 2H for n 1 and provide an explicit expression for Ch- 
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INTRODUCTION 



Long-memory stochastic processes are ubiquitous in fields as different as condensed mat- 
ter, biophysics, social science, climate change, finance The development of 
methods able to quantify the statistical properties and, in particular, to extract the Hurst 
exponent of long-range correlated signals continue therefore to draw the attention not only 



of the physicist community 
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For long-memory corre- 



lated processes, the value of the Hurst exponent H ranges from < H < 0.5 and from 
0.5 < H < 1 for negative and positive persistence respectively; H = 0.5 is found in fully 
uncorrelated signals. Several techniques have been proposed in the literature to study the 
scaling properties of time series. We limit ourselves to mention here only a few of them 
such as the seminal work by Hurst on rescaled range statistical analysis (R/S), the modified 
R/S analysis, the multi-affme analysis, the detrended fluctuation analysis (DFA), the peri- 
odogram regression (GPH) method, the (m,k)-Zipf method, the detrended moving average 
analysis (DMA). The challenge is to get the Hurst exponent H, that is related to the fractal 
dimension D = 2 — H, by means of more and more accurate and fast algorithms. The meth- 
ods of extraction of the scaling exponents from a random signal exploit suitable statistical 
functions of the series itself. The rescaled range statistical analysis (R/S) was first intro- 
duced by Hurst to describe the long-term dependence of water levels in rivers and reservoirs. 
It provides a sensitive method for revealing long-run correlations in random processes. This 
analysis can distinguish random time series from correlated time series and gives a measure 
of a signal roughness. The Detrended Fluctuation Analysis (DFA) is a very popular scaling 
technique to estimate the power-law correlation exponents of random signals in the time 
domain. Recently, a method called Detrended Moving Average (DMA) technique for the 
analysis of the persistence has been proposed [13]. The striking difference between the DMA 
and the previously proposed (R/S, DFA) methods is that a division of the series in boxes is 
not needed by the DMA algorithm. The scaling property is obtained by using a simple con- 
tinuous function: the moving average. This fact makes the DMA algorithm highly efficient 
from the computational point of view. Many properties of the DMA have been studied and 
applications have been demonstrated jl4|. 

In this work, we provide an analytical derivation that the DMA variance scales as a power 
law with exponent 2H, i.e. cr 2 DMA ~ Cnn 2H . The analytical findings are compared with the 
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simulation data obtained by calculating the DMA variance of surrogate Fractional Brownian 
motions with different values of H, generated by the Random Midpoint Displacement (RMD) 
algorithm r|. 

METHOD 

First we describe the main steps of the DMA algorithm. The technique is based on the 
function: 

1 N 2 

i=n 



1 n 



Vn\V 
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The Eq. (jlajl defines a generalized variance of the random path y(i) with respect to the 
moving average y n (i) (Eq. (jlbj)) . The function y n (i) is calculated by averaging the n-past 
value in each sliding window of length n. In so doing, the reference point of the averaging 
process is the last point of the window. The dynamic averaging process and the DMA 
algorithm can be however referred to any point within the window, by generalizing the 
Eqs. (Ilallbj) as follows: 

^ N-nd 2 
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Upon variation of the parameter 9 in the range [1, 0], the reference point of y n (i) is accord- 
ingly set within the moving window n. In particular, consider the following three relevant 
cases: (i) 9 = corresponds to calculate y n (i) over all the past points within the window 
n; (ii) 9 = ~ corresponds to calculate y n (i) over n/2 past and n/2 future points within the 
window n; (iii) 9 = 1 corresponds to calculate y n (i) over all the future points within the 
window n respectively. 

In order to calculate the Hurst exponent of the series, the DMA algorithm is implemented 
as follows. The moving average y n (i) is calculated for different values of the window n, with 
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n ranging from 2 to a maximum value n max determined by the size of the series. The a dm A, 
defined by the Eqs. fllal [Tcjl . is then calculated for all the windows n over the interval [n, N]. 
For each y n (i), the value of <Jdma corresponding to each y n {i) is plotted as a function of n 
on log-log axes. The most remarkable property of the log-log plot is to exhibit a power-law 
dependence on n, i.e. & 2 DMA ~ Cnn 2H , allowing thus to calculate the scaling exponent H of 
the signal y(i). 

DERIVATION OF THE SCALING RELATIONSHIP 

The variance defined by the Eqs. (0) will be analytically deduced in the limit iV 3> n 3> 1 
for fractional Brownian path. The calculation will be performed for the DMA variance with 
moving average y n {i) referred to an arbitrary point within the moving window. We will 
prove that: 

°Ima ~ C H n 2H , n » 1 (2) 

with 
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By simple transformation of the Eq. (fTcj) . one obtains: 
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Let us consider each term separately. The first term on the right hand side of the equation 
(HJ) writes: 

N-9n N~9n 
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The second term on the right hand side: 
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+ ^- 1 [(l-9) 2H+1 + 9 2H+1 ](N-n) (6) 



Finally, the third term: 
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Summing the contributions from each term we obtain 
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One can easily check that <J 2 DMA takes the following expressions for 6 = and 6 = 1 and for 
6 = 1/2: 

• 6 = or 6 = 1; i.e. when the moving average is referred to the last or to the first 
point of the window, it is 



6 = | : the average is referred to the center of the window 

r,rr ~ [ 

■2(H 



1 I _ o~ 2H 

c ^[wttt-wtt] < 10 > 



The above calculations have been performed for a fractional brownian motion with variance 
a = t 2H . It can be easily shown that for the general case of a fractional brownian motion 
with variance a 2 = Dnt 2H , the Eq. (JHJ) will asymptotically behave as cr 2 DMA = DhChit 211 . 
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Intercept (A) and slope (B) for the curves 


plotted in Fig.l. 
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RESULTS AND DISCUSSION 

In this section, the values of Ch obtained by calculating the DMA variance of artificial 
fractional Brownian motions are compared with those calculated using the Eq. (JSJ). 

In Fig. (£[]), the results of the DMA algorithm implemented over artificial fractional random 
walks generated by the random midpoint displacement algorithm are shown. The length of 
the walks is = 2 23 and the Hurst exponent ranges from 0.1 to 0.9 with step 0.1. The 
curves in the three figures refer respectively to three values of the parameter 6, namely 
9 = 0, 9 = 0.5 and 9=1. The slopes of the curves plotted in Fig.((TJ) are plotted in Fig.(j2J). 
The slopes and the intercepts of the curves plotted in Fig.(P) are reported in Table I. From 
the data shown in Table I, it is possible to deduce that the DMA with 9 = 0.5 performs 
better with positively correlated signals with 0.5 < H < 1, while the DMA with 9 = and 
9=1 performs better with negatively correlated signals with < H < 0.5. In Fig.Q, the 
theoretical values of Ch, calculated by using the Eq.©, are compared with those obtained 
by the intercepts of the curves plotted in Fig.(l) for 9 = 0.5 (data of the 2nd column of 
Table I Q). 

It is interesting a comparison of the Eqs. ()9|10j) with the corresponding ones obtained 
for the DFA method, that we will briefly recall here. According to the DFA method, the 
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FIG. 1: Log-log plot of the function a dm A denned by the Eq. ^for artificial series generated by 
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FIG. 2: Values of the slopes of the function ctdma plotted in Fig.Q for three different values of 
the parameter 6 respectively equal to 0, 0.5 and 1. 

integrated profile y(i) is divided into boxes of equal length n. In each box, the signal y(i) 
is best-fitted by an £-order polynomial y n ,e{i), which represents the local trend in that box. 
The different order of the DFA-£ (e.g.,DFA-0 if £ = 0, DFA-1 if £ = 1, DFA-2 if £ = 2, etc) is 
obtained according to the order of the polynomial fit. Finally, for each box n, the variance: 

1 - 

i=l 

is calculated. The calculation is then repeated for different box lengths n, yielding the 
behavior of udfa over a broad range of scales. For scale-invariant signals with power-law 
correlations, the following relationship between the function odfa and the scale n holds: 

°l FA ~n 2H . (12) 



DFA-0 and DFA-1 functions has been investigated by 
jis| and by Taqqu, Teverovsky and Willinger in [l2j 



The asymptotic behavior of the DFA-0 and DFA-1 functions has been investigated bs 
Raymond and Bassingthwaighte in 



respectively. Raymond and Baasingthwaighte worked out the Mowing relation (E q .(21) of 
Ref. |18l |) for the asymptotic scaling of the (Jdfa-o function [20]: 



^ ( k 2F+T ~ 2(F+i) J n8H (13) 
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FIG. 3: Values of Cjj for the DMA algorithm for 8 = 0.5 (squares) and for the DFA-1 algorithm 
applied to the same series (circles). The solid lines represent the values of Ch calculated by using 
the expressions (jHJ and (|14|) respectively. 



The asymptotic behavior of the DFA-1 function is reported in the Appendix of 
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2H + 1 H + 2 H+l / 

In Fig. |3J the values of Ch for the DMA and the DFA are shown. It can be observed that 
the behavior of Ch obtained from the simulations (square and circles) follows quite well the 
analytical curves (solid lines) around H ~ 0.5. Deviations are observed at the extrema of 
the H range. Such deviations might be related either to the DFA and DMA sensitivity or 
to the RMD accuracy. This issue, which is beyond the scope of the present effort, will be 
addressed in a future paper. 

CONCLUSIONS 

We have analytically derived the asymptotic scaling behavior of the DMA algorithm 
proposed in We have compared the analytical results with those obtained from the 

simulations of fractional Brownian paths with assigned values of H obtained by the Random 
Midpoint Displacement. 
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of the RMD algorithm. 

It „ t o cll e ck tb. ^ widowed vaHa.ee w it W a ny treIld in Ref.Q 

is indeed equivalent to the DFA-0 variance. Note that the sentence "without trend correction" 
is not fully appropriate. The function obtained by fitting the random walks y{i) by constant 
segment in each box corresponds indeed to a zeroth-order approximtion of the trend of y(i). 
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